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We propose a scheme for state-insensitive trapping of neutral atoms by using light with two inde- 
pendent wavelengths. In particular, we describe the use of trapping and control lasers to minimize 
the variance of the potential experienced by a trapped Rb atom in ground and excited states. We 
present calculated values of wavelength pairs for which the 5s and 5p3/2 levels have the same ac 
Stark shifts in the presence of two laser fields. 
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I. INTRODUCTION 

The emerging advantages of using transitions between 
different atomic electronic configurations for frequency 
standards and quantum information processing are ac- 
companied by disadvantages with respect to previous 
methods in which atomic qubits were hyperfine states 
of the same configuration. In the latter case, the energy 
shifts of states induced by external trapping fields are 
small, and can be calculated to an accuracy that usually 
does not make a significant contribution to the uncer- 
tainty budget. When qubits are associated with different 
configurations, on the other hand, the magnitude and 
even the sign of differential field shifts are uncontrolled 
in the first instance. 

The ability to trap neutral atoms inside high-Q cavities 
in the strong coupling regime is of particular importance 
for quantum computation and quantum communication 
schemes, where it is essential to precisely localize and 
control neutral atoms with minimum decoherence. In a 
far-detuned optical dipole trap, the potential experienced 
by an atom in its ground state can be either attractive 
or repulsive, with respect to the location of peak light 
intensity, depending on the sign of the ac Stark shift due 
to the trapping light. Excited-state atoms in the same 
trap may experience an ac Stark shift with an opposite 
sign, which affects the fidelity of experiments in which 
excited states are temporarily occupied, such as the im- 
plementation of the Rydberg quantum gate [l|-Q . 

The same problem, i.e. different Stark shifts of 
two states, affects optical frequency standards based on 
atoms trapped in optical lattices, because it can intro- 
duce a significant dependence of the measured frequency 
of the clock_ transition upon the lattice wavelength. Ka- 
tori et al. [^ proposed the idea of using a trapping laser 



tuned to a magic wavelength, A 



magic 1 



at which the ac 



Stark shift of the clock transition is eliminated. The 
magic wavelength of the ^^Sr ^Sq— '^Pq clock transition 
was found to be 813.5±0.9 nm in Ref. [7] by investigating 
the wavelength dependence of the carrier linewidth. This 
magic wavelength was later determined with even higher 
precision to be 813.42735(40) nm [8,]. In a cavity quan- 
tum electrodynamics experiment, McKeever et al. [9| 
demonstrated state-insensitive trapping of Cs atoms at 
A 



magic 
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935 nm while still maintaining a strong cou- 
pling for the 6p3/2 — 6si/2 transition. 

Magic wavelengths for np — ns transitions in alkali- 
metal atoms from Na to Cs have been previously calcu- 
lated by Arora et al. [lu|, using a relativistic all-order 
method. This was accomplished by matching the ac po- 
larizabilities of the atomic npj and n.3 states. The data 
in Ref. |l(j] provide a wide range of magic wavelengths for 
alkali-metal atoms. In the case of the np^/2 ~ i^s transi- 
tions, the magic wavelengths need to be determined sep- 
arately for the rrij = ±1/2 and rrij = ±3/2 states, due 
to the rank-2 tensor contribution to the polarizability 
of the np^/2 level. Furthermore, there is a substantial 
reduction in the number of magic wavelengths for the 
nij — ±3/2 states due to selection rules for linear polar- 
ization. For instance, three out of the six values of Amagic 
suggested for the 5^3/2 — 5s transition in Rb are present 
only for the rrij ~ ±1/2 states. In such cases, the magic 
wavelength becomes dependent on the particular hyper- 
fine state of the atom. Some of the magic wavelengths 
are also in regions that are inconvenient for present laser 
technology. Out of the remaining three wavelengths con- 
sidered in [lO[, the Amagic at 791 nm has opposite signs 
for the Stark shifts for rrij = ±3/2 and rrij = ±1/2 
states, which makes this wavelength of limited practi- 
cal use. The second magic wavelength at 776 nm is in 
close proximity to the Rb 5p — 5d resonance transition 
at 775.8 nm, which could mediate undesired two-photon 
transitions. The third magic wavelength at 637 nm exists 
for all states, but its corresponding polarizability is too 
small for convenient trapping (it is -500 Cq, where uq is 
the Bohr radius). In summary, the single-laser scheme 
offers few cases in which the magic wavelengths are con- 



venient for state-insensitive trapping of Rb atoms [l3| ■ 

In this paper, we investigate an obvious mechanism 
for remediating uncontrolled frequency shifts in transi- 
tions between different configurations: the application of 
a second, "control," optical field to a system of optically- 
trapped atoms. We outline the general principles of this 
approach, and apply it in detail to some cases in Rb. Ru- 
bidium is chosen because it offers a baseline of compari- 
son with previous, monochromatic, attempts at control, 
and this serves to illustrate advantages of the bichromatic 
approach that we believe will have wide applicability. 

Specifically, we find the combinations of two wave- 
lengths that allow to match ac Stark shifts of the Rb atom 
in 5s and 5p3/2 states. In this scheme, a combination of 
trapping and control lasers allows one to minimize the 
difference in the trapping potentials experienced by the 
atom in ground and excited states. This approach sig- 
nificantly increases the number of wavelengths at which 
state-insensitive trapping experiments can be conducted. 
The first step in the realization of this scheme is to 
calculate the Stark shifts of the 5s and hp^/2 states of 
the Rb atom as a function of frequency. We use the 
relativistic all-order method |1CHL3 | for the calculation 
of reduced electric-dipole matrix elements involved in the 
evaluation of frequency-dependent polarizabilities. In the 
second step, we calculate the shift in energies of atomic 
states as a function of the two laser frequencies. The 
wavelengths are determined where the ac Stark shifts of 
the 5s and 5^3/2 levels match according to criteria that 
are described below. Several specific cases are illustrated 
in detail. 



where C is given by 



C = 



^jvC^jv - 1) 



6(j.-Hl)(2j„ + l)(2j„ + 3) 



1/2 



(3) 



and (fc||d||?;) are the reduced electric-dipole matrix el- 
ements. The experimental energies Ei of the dominant 
states contributing to this sum have been compiled for 
the alkali atoms in Refs. [Ij, |T5|- In addition to the 
scalar and tensor valence contributions, there is a scalar 
core contribution to the polarizability, acoro- For the fre- 
quency range considered in this work, Ofcorc has a very 
small uj dependence. The static core polarizability value 
calculated using a random-phase approximation [16[ has 
been used in our calculations without loss of accuracy, 
i.e. uncertainty of this term gives negligible contribution 
to the total uncertainty. 

Unless stated otherwise, we use atomic units (a.u.) for 
all matrix elements and polarizabilities throughout this 
paper: the numerical values of the elementary charge, 
e, the reduced Planck constant, h — h/2'K, and the 
electron mass. Trie, are set equal to 1. The atomic 
unit for polarizability can be converted to SI units via 
a/h [Hz/(V/m)2]===2.48832xl0~^a (a.u.), where the con- 
version coefficient is 47reoag//i and the Planck constant 
h is factored out in order to provide direct conversion 
into frequency units; ao is the Bohr radius and cq is the 
electric constant. 
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II. FREQUENCY-DEPENDENT 
POLARIZABILITY 

The second-order energy shift /S.E of a monovalent 
atom in a state v is parameterized as the sum of scalar 
ao(w) and tensor 02(1^) polarizabilities 



AE = — anfwje a2[uj) —, ^ 



(1) 



where the laser frequency uj is assumed to be several 
linewidths off-resonance, jy is the angular momentum, 
e is the rms magnitude of the electric field, and the po- 
larization vector of the linearly polarized light defines 
the z direction. The valence contribution to frequency- 
dependent scalar and tensor polarizability is evaluated 
as the sum over intermediate k states allowed by the 
electric-dipole transition rules [1^ 
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FIG. 1: Surface plot for the 5s and 6^3/2 rrij — ±1/2 state 
polarizabilities as a function of laser wavelengths Ai and A2 
for equal intensities of both lasers. 

The ground and excited state ac polarizabilities of 
the alkali-metal atoms were previously calculated accu- 
rately in Refs. [13, [13, [13, [ll • Detailed description of 
the polarizability calculations for atomic Rb is given in 
Refs. [10|,[3- Briefly, the sums over intermediate states k 
in the formulas are separated into a dominant part amain 
that contains the first few terms and a remainder ataii- 
In our Rb calculations, we include all ns states up to 10s 



TABLE I: Magic combinations of the trap and control wavelengths Ai, A2 = 2Ai for the 5p3/2mj — 5s transition in Rb and the 
corresponding sum of polarizabilities at these wavelengths. The wavelengths (in vacuum) are given in nm and polarizabilities 



are given in atomic units, e^ and £2 represent the intensities of the two laser beams, respectively, a^ 
that the energy level shift is proportional to a'^^'^ei. 



ai + (e2/ei) a2, so 
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FIG. 2: Surface plot for the 5s and 5p3/2 rrij = ±3/2 state 
polarizabihties as a function laser wavelengths Ai and A2 for 
equal intensities of both lasers. 



set consists of 70 splines of order 11 for each value of the 
relativistic angular quantum number k. 

In the calculation of the main term, the 5p3/2 ~ 5s ma- 
trix elements are taken from Ref. 12 C| , and the 5p3/2 — 4:dj 
El matrix elements are taken to be the recommended val- 
ues derived in Ref. [21| from the Stark shift measurements 
reported in [22] . We use the all-order method (linearized 
version of the coupled cluster approach), which sums in- 
finite sets of many-body perturbation theory terms, for 
the calculation of all other matrix elements in the dom- 
inant part, amain- Detailed description of the all-order 
method is given in Refs. [ll|, [2^. For some matrix ele- 
ments, it was possible to carry out semi-empirical scal- 
ing of the all-order values to include some additional 
important higher-order corrections. The scaling proce- 
dure has been described in Refs. [11,, 13, 24]. The result- 
ing frequency-dependent polarizabilities are used to find 
convenient combinations of trap and control laser wave- 
lengths that yield the same ac Stark shift for Rb atoms 
in the ground and excited 6^3/2 levels. 



and all nd states up to 9d in the amain term. The ataii 
contribution is calculated in the Dirac-Fock (DF) approx- 
imation. We use a complete set of DF wave functions on a 
nonlinear grid generated using B-splines [19| constrained 
to a spherical cavity. A cavity radius of 220 oq is chosen 
to accommodate all valence orbitals of amain- The basis 



III. RESULTS 



In this section, we list a few appropriate combinations 
found for control and trap laser wavelengths where the 5s 
and 5p3/2 state polarizabilities of Rb are closely matched. 



For monochromatic light, a magic wavelength is repre- 
sented by the point at which two curves, a5s(aj) and 
a5p(ci;), intersect as a function of the frequency, u. In 
the bichromatic case, on the other hand, we have two 
additional degrees of freedom, the control frequency and 
the ratio of laser intensities. Thus, bichromatic magic 
wavelengths are represented as curves resulting from the 
intersection of surfaces. 

To illustrate this point, we display sample cases of 
such surface plots for the ttlj — ±1/2 and rrij — ±3/2 
states in Figs. [T] and [2 respectively. The intensity of 
both lasers is taken to be same, so that the energy level 
shift is proportional to the sum of two polarizabilities. 
This is plotted on the z axis. The trap and control laser 
wavelengths are given on the x and y axes, respectively. 
The total polarizability of the 6^3/2 state depends upon 
its rrij quantum number, and it is calculated as a sum 
or difference of the scalar a^ and tensor a2 polarizabili- 
ties, i.e. 0(5^3/2) = ao — q;2 for jtij = ±1/2 states and 
0^(5^3/2) = oiQ + 0L2 for nij = ±3/2 states. Therefore, 
we discuss the results for rrij = ±1/2 and rrij = ±3/2 
states separately. We also find some appropriate trap and 
control laser wavelength combinations that have similar 
magic wavelengths for each rrij state. As illustrated by 
Figs. [Hand [H there is a large number of possible combi- 
nations of trap and control wavelengths that will result 
in the same ac Stark shift of both levels. 

In Table m we list a number of sample trap and control 
wavelength combinations which can be used for state- 
insensitive trapping of Rb atoms in 5s and 5p3/2 states. 
Out of a number of combinations found, we list only those 
where one of the laser wavelengths is twice the other. 
This is a case of particular practical interest, since it is at- 
tainable by frequency doubling of the longer-wavelength 
laser. The combinations are listed for various trap and 
control laser intensity ratios as indicated to illustrate the 
ability to tune the magic wavelength pairs by varying the 
relative intensities. 

The percentage difference between total polarizabilities 
of the 5s and 5ps/2 states for the cases listed in Table H] 
is less than 1 % taking into account uncertainties 

We discuss the magic wavelengths for the nij — ±1/2 
states first. For equal intensities, the combination with 
Ai = 788 nm and A2 = 1576 nm, may be particu- 
larly useful since the resulting polarizability is positive 
and the atoms in red detuned traps are attracted to- 
wards the maximum of the field intensity [25l l26j . Ap- 
plying a control laser with double the trap laser wave- 
length creates a deeper trapping potential for the atom 
in the ground state and minimizes the difference between 
the Stark shifts for ground and excited states. For the 



5P3/2 



1/2 state the polarizability is negative at 



788nm (—10279 Oq) and larger and positive at 1576 nm 
(15194 Oq). The uncertainties in the polarizabilities at 
combinations which are close to resonance wavelengths 
are generally higher. The value of a''"™ for some of the 
combinations in Table U is negative, so that the atoms 
become low-field seekers. A number of groups have sug- 
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3: Magic wavelength pairs for Ai = 800 — 830 nm and 
1600 — 1660 nm and equal intensities of both lasers. 



gested blue detuned or dark optical traps where atoms 
are surrounded by repulsive light fields and; therefore, 
are captured in dark regions without light [2y, [23j . In 
contrast to the monochromatic case, where very few con- 
venient magic wavelengths were found for rrtj ~ ±3/2 
states, a number of "dark" magic wavelengths for Rb are 
found in the present bichromatic treatment. 

We also found a few laser wavelength combinations 
that support state-insensitive simultaneous trapping for 
all iTij states. Examples of such cases are grouped to- 
gether in a first few rows of Table HI The magic wave- 
length combination for \mj\ = 1/2 case is given first, and 
the corresponding |r7T,j| — 3/2 magic wavelength combi- 
nation is given in the following row. We illustrate the 
example of such magic wavelength combinations (listed 
in rows 3 and 4 of Table |l| in Fig. [3l where we plot sur- 
faces of the 5s and 5p^/2 \mj\ = 1/2,3/2 state polariz- 
abilities for Ai — 806 — 826 nm and A2 — 1615 — 1645 nm. 
The intensities of both lasers are taken to be equal. The 
magic wavelength by may be further tuned by adjusting 
intensity ratio of the two lasers as illustrated in Fig. S) 



The magic wavelength for the 5s and 5p3/2 



1/2 



states for Ai = 800 — 810 nm and A2 = 2Ai are shown 
for various intensities of both lasers. The intensity ratio 
(£1/62)^ ranges from 1 to 2. 



IV. CONCLUSION 

In summary, we have explored a bichromatic scheme 
for state-insensitive optical trapping of Rb atom. Due 
to the extensive development of first principles atomic 
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structure theory, semiempirical corrections, and compu- 
tational methodology we are able to explore a wide range 
of parameter space with reasonable confidence in the un- 
certainties of our calculations. We have recently com- 
pleted a comprehensive survey of calculations of DC po- 
larizabilities, for which there exist copious experimental 
data for comparison within clearly-stated ranges of un- 
certainty [28|. In this paper, we specifically explored a 
case of the Rb atom, where the magic wavelengths as- 
sociated with monochromatic trapping were sparse and 
relatively inconvenient. We have found that the bichro- 
matic approach yields a number of promising wavelength 
pairs which are discovered with straightforward parame- 
ter choices such as equal laser intensities and A2 = 2Ai. 
The methodology developed in this work allows us to ex- 
plore specific cases of interest that may arise in the future 
experiments where it is essential to precisely localize and 
control neutral atoms with minimum decoherence. 



FIG. 4: Magic wavelength for the 5s and 5p3/2 rrij = ±1/2 
states for Ai = 800 — 810 nm and A2 = 2Ai for various inten- 
sities of both lasers. The intensity ratio (ei/e2)^ ranges from 
1 to 2. 
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